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Abstract 

A method is shown to estimate the position of by fusing 
the data from a sodium iodide detector and a gyroscope while 
panning the detector. Based on simple geometry considerations, 
the search motion of a detector is modelled as angular panning. 
Correlation of both sensor streams is shown as a predictor for 
the direction of the incoming radiation. The method also allows 
for a rough distance classification, where a distant source or a 
homogenous radiation field can be distinguished from a single 
source near to the detector. 


1 Introduction 

Handheld radioisotope identification devices (RIID) are widely used as security 
equipment. They complement static portal detection philosophies and allow 
their users to actively search for sources by moving the device. Most instruments 
therefore contain an extra mode for this search to give visual, acoustic and 
vibrational feedback about the strength of the current radiation field. In the 
best case, they finally guide their users to the hotspot or source of radiation. 

A very interesting feature would be a compass-like needle, pointing roughly 
towards the direction where a radiation hotspot is assumed. Please regard also 
to [I] and reference therein. There are different ways to do this. First, one 
can use multiple detectors, in the easiest configuration two. Then, common 
approaches as described in [2] apply and the position of the source can be 
approximated even with a certain amount of perturbation material within the 
path of the incoming radiation. But there are also several reasons not to use 
more than one detector in an handheld instrument: (a) it is expensive, (b) it 
increases the complexity of the electronics system and (c) it increases the points 
of possible failures. 

More elaborated approaches also use multi-detector concepts such as coded- 
aperture applications 0111,0 or even Compton-cameras. The latter have seen 
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a significant evolution and with statistical algorithms as presented in [6], such 
cameras can locate a source very reliable. Nevertheless, solutions like these two 
examples are not within our scope. We refer to a rather inexpensive handheld 
configuration with only one detector system. 

We therefore adopted a different approach based on sensor or data fusion 
[7], whereas we explicitly refer to the merging of information from two different, 
heterogenous sources. In our case this is the nuclear detection system and a 
gyroscope (gyro) [8]. Adding a gyro to the detection system is an easy task. 
Such modules are custom-out-of-the-shelf products and can be coupled to a lot 
of interesting hardware concepts, e.g. a Raspberry Pi micro computer platform. 
On the other hand, most mobile platforms like Android boards or the iPhone 
are already fully equipped with these sensors. Once the data of the gyro is 
analysed, we can reconstruct a part of the local detector motion. Together with 
the count rate given by the radiation detector we have enough data to predict 
the direction to the source. 

Data fusion has several prospects as it adds information for the interpretation 
of a given sensor. Nowadays, it is used in a wide range of applications [9j 
Eg. Most prominently though not stated exactly as data fusion, examples also 
include the application of sensor fault detection m ■ 

Our technique represents a virtual or so-called soft sensor. With that, the 
fused result of both sensorial inputs represents a new, algorithmically deduced 
sensor component that gives the answer to the questions from where the source 
radiation is coming from. 

Note that we do not involve any shielding of the detector, nor do we restrict 
to collimated source beams. Instead we want to establish a simple, robust but 
certainly rough mean to estimate the direction to the source. 

Our presentation starts with some very basic models to gain insight into 
the setup. It is clear to us, that real field missions are practically unforeseeable. 
The multitude of different measurement geometries not only including scattering 
and absorption effects, but also distributed sources, render it quasi impossible 
to calculate the radiation direction precisely. Nonetheless, we believe that a 
rough estimator, e.g. indication whether the source is left or right of the search 
path, is still achievable. 


2 Sensor fusion approach 

2.1 Approximative relation between count rate and geom¬ 
etry 

Key idea is the hereby the knowledge of the gyro states, in detail knowledge of 
the local angular motion of our instruments. We record this motion in terms of 
three coordinates <a,/3,y, their respective velocities v a = da/dt , vp = d/d/dt, v 1 = 
d^/dt and accelerations a a = d 2 a/dt 2 , ap = d 2 /d/dt 2 , a 7 = d 2 ^/dt 2 . Available 
gyro components grant experimental access to these quantities. For the sake of 
simplicity, we define the angle a to correspond to our motion. That means, we 
pan the instrument around the z-axis, within the x - 7/-plane. 

We apply a very simple circular model for the search motion of a handheld. 
Lets for a moment assume we know the distance D between instrument rotation 
center and source. R may be the radius of our search motion, a the angle 
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Figure 1: Sketch of our geometry, (a) Side view of the detector and our axis 
definitions as well as the angle a. (b) Top view. Distance of the rotation center 
to the source D. The detector panning is assumed to be circular. 


reported by the gyroscope, with d = D - Rcosa it follows that 

Ar 2 = (D - Rcosa) 2 + R 2 sin 2 a. (1) 

Here x is now the resulting distance to the source during rotation. Please 
regard to Fig. [l]for further details. 


2.2 Count rate 

In our model of the source we define the count rate at the sources physical surface 
to be Ao. Following the quadratic law, we can relate the distance between source 
and detector to the count rate [I], 


A(Ar) 


Ap 

(Ar + l) 2 


and further to angle a from 0, 


A(a) 


_ Ap _ 

((D - Rcosa) 2 + R 2 sin 2 a + l) 2 


( 2 ) 


(3) 


2.3 Simplified models of detector motion 

2.3.1 Rotation 

We first start with a simple rotation. It can be easily checked whether the 
output of our gyroscope considerations is correct, once the In Fig. [2] we show an 
analytical evaluation of yielding our assumption on what the count rate at 
our detector should look like, if it is rotated along the z -axis within the radiation 
field. The result depends on the distance D and the panning radius R. Here it 
is natural that with D/R -> oo the impact of rotating completely vanishes. 
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Figure 2: Count rate A as function of the angle a , while the detector would be 
rotating. Different ratios D/R are plotted. 


2.3.2 Waving 

The most relevant search motion for finding a source is panning or waving, in 
other words moving the detector from right to left and back again. With the 
previous Eq. ©> this behaviour can be described by restricting the circle to sort 
of waving cone. In right part of Fig. [3] we show a sketch of this motion. 

Following scenarios can be found with this motion and for each scenario 
Fig. 1 depicts the expected count rate as function of time: (a) the source is at 
the left-(right)-hand-side of the detector. In this situation the count rate rises 
till we reach the point closest to the source. It would be a maximum, if the 
source is really in line with the final position of the detector, (b) the source 
within the 90° waving cone, yielding an elongated maximum in the count rate, 
(c) the source is approximately at the middle of the left quarter of our waving 
motion. The count rate rises, has a first maximum when passing the source 
direction, decreases to an intermediate minimum, because the source is not as 
far away as on the other side and finally it passes the source direction again 
with the returning motion leading into the next oscillation. This curve and the 
associated source position are marked in red. As last case, (d) marks a situation 
where the incoming radiation is close to the center axis of our waving motion 
and no indication is possible. Here, the count rate yields an oscillation peaks 
with similar heights. Note that here, each panning motion passes the highest 
count rate twice, so there are also twice the number of peaks as in scenario (a). 

2.4 Forecasting based on the angular information 

Of course, the ratio D/R has the same impact on the results as in the case 
with rotation. From the shapes found by panning the detector we constitute 
following data fusion rules: 

• Once there is no change of count rate during panning, there are distributed 
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Figure 3: Count rate A as function of the angle <a, while the detector is panned 
for the ratio D = 2 R and an approximate count rate of A = 1600cps. For compar¬ 
ison with real sources, our detector has a sensitivity of about S = 2400cps//iSv/h 
for cesium 137 Cs. The simulation assumes a distance so that the dose rate be¬ 
comes p = 1.6/iSv/h. 


sources or a single source far away. 

• If the function of count rate yields significant changes, especially peaks, 
there is at least one source close by. 

• The number of oscillation peaks is a first indication whether the incoming 
radiation is coming from the front or from left or right. 

and last, but perhaps our most important statement, 

• Let a(t) be the current local waving angle as function of the time t and let 
a(t) be the count rate of our detector, then the correlation value between 
both C\a(ts)i is an directional measure for the incoming radiation. 


3 Implementation in the prototypical and com¬ 
mercial scope 

3.1 Prototype implementation 

A simple prototype was assembled using a Raspberry Pi micro-computer and 
an I 2 C based gyro module. A multi-channel analyser was mounted on a 3" x l" 
sodium iodide detector and connected to the Raspberry Pi via serial communi¬ 
cation. To verify our ruleset presented in |2.4| we programmed a python software 
running on the Raspberry Pi, reading out the gyroscope and the MCA counts 
simultaneously. Using the prototype, the statements of |2.4| were confirmed and 
the algorithm was further improved with regard to performance. 
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Figure 4: Count rate A as function of the angle <a, while the detector is panned 
for the ratio D = 2 R and an approximate count rate of A = 1600cps. For 
comparison with real sources, our detector has a sensitivity of S = 980cps//iSv/h 
at the energy of the cesium 137 Cs photo peak at 661.65keV. The simulation 
assumes a distance so that the dose rate becomes p = 1.6/iSv/h. 

3.2 Commercial exploitation 

To demonstrate the usability and robustness of the developed technique, we in¬ 
tegrated a gyroscope into the hardware of a commercially available radioisotope 
identification device (RIID), which features a specific software user interface 
called ” detect mode”. This significantly extended the scope of our initial proto¬ 
type. Figure [4] shows the algorithm behind the directionality. Depending on the 
determined value of C, a visual indication is given to the user of the instrument, 
in which direction the source is located. The instrument automatically identifies 
the motion type by its time pattern. 
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